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1 Introduction 


In a previous paper ||^ we presented a model to explain the origin of the Stan¬ 
dard Model (SM) fermion mass and mixing hierarchy, based on the so-called 
anti-grand unified extension of the SM gauge group (SMG). This model gives 
rise to a characteristic structure for the quark-lepton mass matrices; in this pa¬ 
per we shall try to argue that some features of this model are implied by the 
experimentally determined fermion masses. We shall also extend the analysis in 

to include CP violation due to a complex phase in the CKM quark mixing 
matrix and to calculate the neutrino mass matrix to get predictions for neutrino 
masses and mixing angles. 

In section we will describe how to calculate fermion masses and mixing an¬ 
gles when viewing the SM as an effective theory. The observed fermion masses 
and mixing angles will be discussed in section |^. An approximate parameteri- 
sation of the masses will then be described. In section |4| we will give a simple 
method by which a model could naturally explain such a parameterisation. Our 
anti-grand unified model will be described in section |^, where we show that it 
fulfils the requirements of section ^ The gauge quantum numbers of the quarks 
and leptons in our model are fixed by the requirement of anomaly cancellation. 

In section ^ we will describe the construction of our model for the second 
and third generation fermion masses and mixing angle, in terms of two Higgs 
fields in addition to the usual Weinberg-Salam (WS) Higgs field. We emphasise 
in section ^ the general features, suggested by phenomenology and the anti¬ 
grand unified model, which underlie the structure of our mass matrices for the 
second and third generation fermions. The choice of the Higgs field quantum 
numbers and the extension of the mass matrices to include the first generation 
particles are discussed in section |[ This leads to definite predictions for the 
order of magnitude of the fermion masses and mixing angles in terms of the 
vacuum expectation values (VEVs) of the Higgs fields. In section ^ we present 
our best fit to the conventional experimental masses and mixing angles, in terms 
of just three VEVs of the same order of magnitude. We also make a fit using 
preliminary values for the light quark masses extracted from lattice QCD. These 
smaller “lattice” values for the light quark masses motivated us to consider, in 
section |^, the possibility of an alternative mass matrix structure 

Then we discuss other predictions of our model. In section we consider 
how much CP violation would be expected from the CKM matrix in our model. 
Finally we turn to the neutrino masses and mixing angles in section |l^. The 
effective Majorana mass matrices for the three light neutrinos in models with a 
hierarchical mass matrix can naturally have quasi-degenerate mass eigenstates 
with maximal mixing . Maximal neutrino mixing provides a candidate expla¬ 
nation for the atmospheric muon neutrino deficit or the solar neutrino problem 
with vacuum oscillations. Our model does have such a structure but, with a 
see-saw mass scale equal to the Planck scale, the predicted neutrino masses are 
too small to give observable effects. As in most other models, it is necessary 
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to introduce an intermediate mass scale in order to obtain observable neutrino 
mixing. This is not so attractive in the anti-grand unification model. However 
we can obtain identical structure for the fermion mass matrices in an anomaly 
free SMG 0 U(l)^ model, in which the fundamental scale is unconstrained and 
could be taken as intermediate between the electroweak and Planck scales. We 
present our conclusions in section |l^. 


2 Modelling fermion masses 

In the SM the fermions (apart from the neutrinos) all get a mass via the Higgs 
mechanism. Before electroweak symmetry breaking there are interactions of the 
form: 


Cmass = QlHu^wsUr + QlHd^wsDr + LlHe^wsEr + h.c. (1) 

where <i>rvs is the WS Higgs field, Ql is the 3 SU(2) doublets of left-handed 
quarks, Hjj is the 3x3 Yukawa matrix for the up-type quarks, etc. If we 
represent the SU(2) doublets <i)rvs and Ql as 2 component column vectors, we 
define: 

= ( 1 0 ) 

and _ 

(3) 

where Ul are the CP conjugates of the 3 left-handed up-type quarks. 

After electroweak symmetry breaking the WS Higgs field gets a vacuum 
expectation value (VEV) and we would write: 

-Cmass = UrMijUr + DrMrDr + ErMrEr + h.C. (4) 


where the mass matrices are related to the Yukawa matrices and WS Higgs VEV 
by: 


M = H 


{4>ws) 

V2 


(5) 


and we have chosen the normalisation so that: 

{cj^ws) = 246 GeV (6) 


The masses of e.g. the 3 up-type fermions are obtained from Mr by diago- 
nalising the matrix to find the 3 eigenvalues. In particular we can find unitary 
matrices Vu and Vd so that: 

V^MuMIjVu = Aia.g{rnl,ml,rnl} 

= diag{m^,m2,TOfc} 
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(7) 

( 8 ) 




The quark mixing matrix is then defined as: 


VcKM = V^Vd 


( 9 ) 


This is where we find aesthetic problems with the SM. First, we have no 
way of calculating the mass matrices since there are no constraints due to e.g. 
gauge symmetries. The elements of the 3 Yukawa matrices are allowed to be 
arbitrary complex numbers. This means that there is no understanding of the 
origin of the masses or mixing angles within the SM. The second problem is that 
we would expect that, since there is no distinction between the 3 generations, 
the order of magnitude of the mass of each fermion would not depend on which 
generation it was in. Further, since, as far as the Higgs sector is concerned, there 
is no distinction between the different types of fermions, we would generally 
expect that all fermions in the SM (except for the neutrinos) would have the 
same order of magnitude masses. Indeed naturality would suggest that all the 
Yukawa matrix elements be of order unity and the fermion masses of order 
{(l>ws)/V^ = 174 GeV. This is clearly not the case. 

So in order to understand the masses of the fermions we must postulate 
some model beyond the SM, which can give different SM fermions very different 
masses without the need for arbitrarily small Yukawa couplings. One method 
to do this is to extend the SM gauge group so that the 3 generations are not 
equivalent under the more fundamental gauge group. Then we would consider 
eq. to be the effective Lagrangian for the Yukawa-Higgs sector. In the full 
gauge group such terms would not generally be gauge invariant since the SM 
fermions would have extra “charges” (generally Abelian and non-Abelian). Of 
course such terms would still be expected in the effective theory at energies 
much below the symmetry breaking scale of the full gauge group to the SMG. 
The difference is that now we have no reason to expect that all the Yukawa 
interactions in the effective theory should be of the same order of magnitude. 

As an example consider the bottom quark. If we suppose that the dominant 
transition between the left- and right-handed components in the fundamental 
theory involves not only the WS Higgs field, but e.g. 2 other Higgs fields W 
and T as shown in fig. |^, we would get the following relation for the effective 
Yukawa coupling of the bottom quark in the SM: 


Mp ]\dp 


( 10 ) 


where (VF) and (T) are simply the VEVs of the new Higgs fields W and T, 
and Mp is the (fundamental) mass scale of the intermediate fermions. Here we 
are assuming that all fundamental Yukawa couplings (the Ai) are of order 1. 
Thus the order of magnitude of the effective SM Yukawa coupling constants are 
given by the product of small symmetry breaking factors, like and So 
now we can explain why the fermions in the SM have different masses and we 
can construct explicit models, by choosing extended gauge quantum numbers 
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Figure 1: Feynman diagram for bottom quark mass in the full theory. The 
crosses indicate the couplings of the Higgs fields to the vacuum. 

for the SM fermions and making specific choices of Higgs fields W, T, etc. to 
allow the transitions in the fundamental theory which give masses to the SM 
fermions. 

It is important to note that this type of model cannot give exact predictions, 
since we are still unable to calculate fundamental Yukawa couplings. Essentially 
we introduce a lot more fundamental Yukawa couplings but then make the 
naturality assumption that they are all of order 1. However, the problem we 
are addressing is the huge range of fermion masses and so we can get some 
understanding of this hierarchy without any knowledge of why the masses are 
exactly what they are. Furthermore we assume that there exists a spectrum of 
vector-like fermion states, all having a mass of order Mp, which can mediate all 
of the required symmetry breaking transitions. 

Clearly there are many different models we could propose to model the 
fermion masses. In order to make progress it is therefore necessary to examine 
the experimentally measured masses and look for relations such as order of 
magnitude degeneracy. One point to note is that the effective Yukawa couplings 
we predict will be the values at the fundamental scale Mp (assumed to be of 
order the Planck scale Mpianck — 10^® GeV) rather than experimental scales 
such as I GeV. 


3 Measured fermion masses 

The masses of the fermions in the SM are usually quoted as running masses, at 
some scale such as I GeV, except for the top mass which is generally quoted 
as a pole mass. These masses can be evolved to other scales by using the 
renormalisation group equations (RGEs). One of the assumptions we make is 
that the SM is valid up to the Planck scale, but assuming a lower fundamental 
scale (say of order lO^® GeV) makes no essential difference. So we evolve the 


5 



Table 1: Yukawa couplings and running masses of the SM fermions at 1 GeV 
and Yukawa couplings at the Planck scale. The top quark pole mass is taken to 
be 180 GeV. 


Fermion 

m(l GeV) 

h{l GeV) 

h{^Mpianck) 


u 

4MeV 

2.3 X 10"^ 

5.4 X 10-3 

-12.1 

c 

1.4 GeV 

8.0 X 10-3 

1.9 X 10-3 

-6.3 

t 

240 GeV 

1.38 

0.39 

-0.9 

d 

9 MeV 

5.2 X 10-3 

1.3 X 10-3 

-11.3 

s 

200 MeV 

1.1 X 10-3 

2.8 X 10-4 

-8.2 

b 

6.3 GeV 

3.6 X 10-2 

7.4 X 10-3 

-4.9 

e 

0.5 MeV 

2.9 X 10-3 

2.9 X 10-3 

-12.8 


105 MeV 

6.0 X 10--^ 

5.9 X 10-4 

-7.4 

T 

1.78 GeV 

1.0 X 10-2 

1.0 X 10-2 

-4.6 


masses to the Planck scale to compare them. The main effect is to change the 
ratios of the quark to lepton masses and also the ratio of the top quark to other 
quark masses. At a scale /r we can express the running masses in terms of the 
Yukawa couplings by: 


m{fi) = h{fi) 


{4>ws) 

V2 


( 11 ) 


The pole mass of a quark is given to leading order in the strong coupling by: 


M = m{M) 




( 12 ) 


In table ^ we show typically quoted values for the SM fermion Yukawa couplings 
at 1 GeV 0 and the corresponding Yukawa couplings at the Planck scale evolved 
using the 1-loop RGBs for the SM (see e.g. [|j). There is some ambiguity, 
particularly for the light quarks, in extracting these Yukawa couplings from 
experiment. We shall consider alternative smaller values for the u, d and s 
quark masses, suggested by recent lattice calculations [^, in section ||. 

In table || we show the magnitudes of the 3 mixing angles (without the GP 
phase) at 1 GeV and the Planck scale. We can see that Vus does not change 
significantly, but the other 2 mixing angles are different at the two scales. We 
find that the running is independent of the CP phase chosen. We also find, in 
agreement with [^, that the ratio of Vcb to Vub is approximately constant. The 
difference at the two scales is only sensitive to the top quark mass. 

If we now look for order of magnitude equalities at the Planck scale we can 
see that, with the exception of the top and probably also the charm quarks, the 
Yukawa couplings of the fermions within each generation are order of magni- 
tudewise degenerate. We could explain the down-type quarks and the leptons 
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Table 2: Experimental values of the mixing angles at 1 GeV and the Planck 
scale. 


Mixing angle 

V{1 GeV) 

Planck) 

ln{V{MPlanck)) 

Vus 

0.22 

0.22 

-1.5 

Vcb 

0.041 

0.049 

-3.0 

Vub 

0.0035 

0.0042 

-5.5 


having degenerate masses in a grand-unified model such as the well known SU(5) 
model. However, in reality, the degeneracy is only true in an order of magnitude 
approximation, whereas we would expect exact equality at some scale in such 
a mode]|^ Another feature, which indicates that the degeneracies are not due 
to unification, is that the up quark is approximately degenerate with the down 
quark and the electron but the charm and, most obviously, the top quarks do 
not fall into the same pattern. Obviously we could say that the up quark is 
the exception and it is simply chance that it has the same order of magnitude 
mass as the down quark and the electron, but this is not entirely satisfactory. 
For example, if the up-type quarks get masses by a different mechanism (e.g. 
supersymmetry with a large tan/3 which could produce the observed top and 
bottom quark masses even if they actually had similar Yukawa couplings), then 
it would seem to be natural to expect that the up quark would be heavier than 
the down quark and the electron, just as the charm and top quarks are heavier 
than the other fermions in their generations. Therefore we take the view that 
there is an approximate mass degeneracy between charged fermions within each 
generation, which is certainly true for the lightest generation. However, in the 
other two generations, the charm and top quarks are exceptions to this rule. 
This means that the charm and top quarks must get their masses by a different 
mechanism from the other fermions. 


4 A Natural Explanation 

If we consider the hypothetical situation where the charged fermions within each 
generation have the same masses, then it would be natural to conclude that the 
masses were generated by some mechanism which didn’t distinguish between 
each type of fermion, but did distinguish between fermions in different gener¬ 
ations. This would indicate that, in some model more fundamental than the 
SM, there should be a distinction between the 3 generations. We shall refer to 
these distinct “generations” in the fundamental theory as “proto-generations”. 
We cannot really claim that the SM fermion masses could be described by an 

^It is of course possible to avoid the predicted degeneracy by introducing a non-minimal 
Higgs structure^. 
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approximation to such a model, since the charm and top quarks are clearly not 
even order of magnitudewise degenerate with the other fermions in their gener¬ 
ations. However, we shall show that the SM fermion masses can be described 
by a very similar type of model. 

An important point to consider is that the fermion masses actually tell us 
only what the eigenvalues of the mass matrices are, they do not specify the 
complete matrices. It seems reasonable to assume that in a model where each 
generation is distinguished (e.g. by having different charges under some new 
gauge interaction), the dominant elements in the mass matrices which would 
lead to the fermion masses could easily be elements on the diagonal, i.e. the 
transitions between left- and right-handed fermions within the same genera¬ 
tion. Now we could explain approximate degeneracy of fermions within each 
generation by a model which produced mass matrices with the same order of 
magnitude diagonal elements. The important point is that we need not require 
the complete matrices to be similar, as long as the off-diagonal elements make 
no significant contributions to the eigenvalues. 

Now we can explain why we have the exceptions of the top and charm quarks. 
To first approximation the largest eigenvalue of a matrix is simply the largest 
element. So in the down-type and charged lepton mass matrices, Mjj and Me, 
the largest element is the one corresponding to the interaction between the 
left- and right-handed fermions in the 3rd proto-generation. In the up-type 
mass matrix, Mjj, it happens that our assumption about diagonal elements 
being dominant was wrong. This is not totally surprising, since we don’t really 
have any good reason to make such an assumption without knowledge of the 
fundamental mass-generating mechanism. Since we know that the left-handed 
top and bottom quarks are in the same SU(2) doublet in the SM, this means 
that the right-handed top quark is not really the fermion in the 3rd proto¬ 
generation of this model. The right-handed top quark must therefore belong 
to the 1st or 2nd proto-generation. This makes no difference in the SM since 
there is no distinction between the generations, but in such extended models 
a distinction does exist. By identifying the right-handed top quark with the 
2nd proto-generation and the right-handed charm quark with the 3rd proto¬ 
generation, we have a simple mechanism to explain why their masses differ from 
the other fermions within their generations. 

So now we have some important criteria for a model of the fermion masses. 
We have made the assumption that the SM is valid up to some high scale, such 
as the Planck scale. We have shown that there are several order of magnitude 
degeneracies when the SM fermion Yukawa couplings are compared at this scale. 
This has led to the requirement that the generations should be distinguished 
in the fundamental model, and that the fermions within each of these proto¬ 
generations should have order of magnitude degenerate masses. However, the 
interactions between generations must distinguish between the charged leptons 
and the 2 different types of quarks, so that the top and charm quark masses can 
be explained as originating from interactions between left- and right-handed 



fermions from different proto-generations (most simply being due to the fact 
that the right-handed top quark of the SM actually belongs to the 2nd proto¬ 
generation and the right-handed charm quark to the 3rd proto-generation). Note 
that this mechanism would not work if, for example, only the top quark mass 
deviated from the “rule” of degeneracy within generations. 

5 The anti-grand unification model 

We will now propose a model which has all the features proposed in the previous 
section. We will first describe the model and then show that it does in fact satisfy 
all our requirements. The model we are considering is called the anti-grand 
unification (AGUT) model. It has previously been considered as a candidate 
for explaining the fermion masses |^. The gauge group for the model is: 

G = SMGi 0 SMG 2 0 SMG 3 0 U(l)/ (13) 

where we have defined: 

SMG, =SU(3)i0SU(2), 0U(1), (14) 

The three SMGi groups will be broken down to their diagonal subgroup, which 
is the gauge group of the SM. The U(l)/ group will be totally broken. So rather 
than unifying the simple factors of the SMG, the AGUT model “splits” them 
into 3 copies of each (hence the name of the model). This gauge group (strictly 
speaking without the U(l)/ group) has been used to successfully predict the 
values of the running gauge coupling constants in the SM at the Planck scale 
[ p^ , as critical couplings estimated using lattice gauge theory. 

We put the SM fermions into this group G in an obvious way. We have 
one generation of fermions coupling to each SMGi in exactly the same way as 
they would couple to the SMG in the SM. The broken chiral gauge quantum 
numbers of the quarks and leptons under the symmetry groups SMGi can readily 
explain the mass differences between fermion generations but cannot explain all 
the mass splittings within each generation, such as the ratio of the top and 
bottom quark masses. It is for this reason that the Abelian flavour group U(l)/ 
is introduced. We then choose U(l)/ charges with the constraint that there 
should be no anomalies and no new mass-protected fermions. This leads us 
almost uniquely to the set of charges yf shown in table ||. We have labelled the 
fermions coupling to SMG^ by the names of the ‘i’th generation of SM fermions. 
However, this is just a method of labelling the representations of the full gauge 
group and, as we have discussed in section we expect that, for example, the 
fermion we have labelled cr will in fact turn out to be the right-handed top 
quark in the SM. 

Although the quantum numbers of the fermion fields are determined by the 
theoretical structure of the model (in particular the requirement of anomaly 
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Table 3: U(l)/ charges of the SM fermions. 


Fermion 

yf 

Fermion 

yf 

Fermion 

yf 

UL 

0 

CL 

0 

th 

0 

Ur 

0 

cr 

1 

tR 

-1 

dR 

0 

Sr 

-1 

bR 

1 

SL 

0 

ML 

0 

tl 

0 

Sr 

0 

Ml: 

-1 

tr 

1 


cancellation), we do have some freedom in the choice of the quantum numbers 
of the Higgs helds. So, to model the fermion masses in the SM, we must choose 
how to break G down to SMG. 

A crucial simplification can be made by considering only the Abelian charges, 
when formulating a model for symmetry breaking. We can justify this by not¬ 
ing that all fermions are in singlet or fundamental representations of the non- 
Abelian groups. In the SM the fermions obey a charge quantisation rule: 

I -h i “duality” -b ^ “triality” = 0 (mod 1) (15) 

where y is the conventional weak hypercharge. Here “duality” is 1 for the fun¬ 
damental representation of SU(2) (the doublet) and 0 for the singlet. Similarly 
“triality” is 1 for the SU(3) triplet, -1 for the anti-triplet and 0 for the singlet. 
Therefore we see that, if these are the only allowed non-Abelian representations, 
we can calculate the non-Abelian representations from the weak hypercharge. 
Similarly in the AGUT model we have such a charge quantisation for the weak 
hypercharge yi associated with each SMGi separately, so that the 4 Abelian 
charges completely specify the representation under the full AGUT gauge group. 
We assume that the Higgs scalar fields also satisfy the charge quantisation rule 
for each SMGi and belong to singlet or fundamental representations of the 
non-Abelian groups. 

As we discussed in section |[ we wish to have a model where the diagonal 
elements in the different mass matrices Mjj, Mo and Me are the same, but 
off-diagonal elements should, in general, differ in different matrices. We can 
now show that this will always be the case in the AGUT model. 

Let us define the U(l) charge vector: 

® = (.T'T'yWj 

Then the net charges of the Higgs helds, other than the WS Higgs held, in a 
transition from one left-handed fermion //, to a right-handed fermion of the 
same type, are given by: 

AQ/g = Qf^ - Qg^ ± Qws (17) 
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where we have a plus sign for up-type quarks and a minus sign for down-type 
quarks and charged leptons. Mass matrix elements with the same value of AQ jg 
will be mediated by the same set of Higgs fields and suppressed by the same 
product of symmetry breaking parameters. To compare elements in the same 
place in different mass matrices we will define: 

AQTij = AQTiTj (18) 

and to consider the diagonal elements we will simplify the notation further by 
defining: 

AQtj = AQtm (19) 

where Ti refers to the fermion in the ‘i’th proto-generation of type T; e.g. D 2 is 
equivalent to s, the strange quark. 

Since we can use the Higgs fields W and (which have opposite charges) 
equivalently in non-supersymmetric models, we will have the same order of 
magnitude diagonal elements if: 

AQc/i = TAQdz = iAQsi (20) 

for any combination of plus and minus signs (possibly dependent on i). From 
eq. (0) we can see that, if we choose both signs to be minus, the charges of the 
WS Higgs field Qws cancel out in the relations of eq. (|2^) . We can then easily 
see from table || that the U(l)/ charges allow: 

AQc/i = —AQoi = —^Qei (21) 

The charges yj for these diagonal elements are all 0 for j ^ i, so the only 
requirement left is that the charges yi satisfy eq. (|^). This is true for the 
simple reason that these are the usual weak hypercharges of the fermions in the 
SM and the relation is just the condition that the fermions have the correct 
charges to couple to the same WS Higgs field. Therefore we have proved that 
the diagonal elements in the 3 mass matrices Mjj, Md and Me are the same 
for each matrix, no matter what Higgs helds we choose or how we extend the 
quantum numbers of the WS Higgs field. Furthermore it is fairly obvious that, in 
general, the off-diagonal elements will not be the same in the different matrices. 
So the AGUT model has satisfied the requirements discussed in section 

6 Constructing a model of the Higgs fields for 
the 2nd and 3rd generations 

There are many ways to approach the construction of a realistic model of the 
fermion masses. In section ^ we argued that there should be approximate degen¬ 
eracy of masses between fermions within each generation, with the exception of 
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the top and charm quarks. We explained how this could be achieved naturally in 
a model where the diagonal elements of each mass matrix were of the same order 
of magnitude but, in general, other elements were different in each matrix. As 
we noted in section |^, this is precisely the case in the AGUT model. Of course 
there are other possible models where the mass matrices are of this form. How¬ 
ever, one important point about the AGUT model is that we have not enforced 
such a condition; rather it was forced upon us by the requirement of anomaly 
cancellation, which determined the charges of the fermions. So, however we 
approach the construction of a definite model (by making a definite choice of 
the Higgs fields to break the AGUT gauge group down to the SMG), there is 
always the possibility of natural order of magnitude degeneracies between the 
masses of different fermions within the same generation. Whether or not this 
actually occurs depends on whether or not the appropriate diagonal elements 
give the dominant contribution to the eigenvalue of the relevant matrices. 

We have previously presented one method for constructing a realistic model 
of the fermion masses and mixing angles |^. Here we will give an alternative 
method which leads to the same model. First we shall choose the quantum num¬ 
bers of the WS Higgs field so that this is the only Higgs field required to produce 
the top quark mass and, thereby, the top quark Yukawa coupling is naturally of 
order one. As we have explained in section |^, in this type of model the dominant 
element in the up-type mass matrix cannot be a diagonal element, since then we 
could not predict that the top quark was much heavier than the bottom quark 
and tau lepton. So we choose the element corresponding to a transition between 
tL and cr, as defined in section ||, to be dominant. Gorrespondingly the element 
which determines the charm quark mass corresponds to the transition between 
Cl and tR. So all we have really done is relabel the right-handed top and charm 
quarks which is simply a matter of definition (identification of fermions in our 
model with those of the SM), since they have the same gauge representations 
in the SM. 

We could in fact choose the dominant element to be any of the 6 off-diagonal 
elements in the up-type mass matrix. However, interchanging SMG 2 and 
SMG3 along with a change of sign of the U(l)/ charges simply corresponds 
to relabelling the second and third proto-generations. This relates the 6 ele¬ 
ments in pairs so that there are only really 3 distinct choices to be made. It 
turns out that, with our assumption that the first generation fermions are all 
order of magnitudewise degenerate in mass and that their masses are due to 
the same diagonal element in the three mass matrices, the 3 choices lead to 
equivalent models for the second and third generation masses. This is because, 
after making the appropriate choice of WS Higgs quantum numbers, the quan¬ 
tum numbers of the transitions among the second and third generations are 
the same in all 3 cases up to permutations of the SMGi and rescaling of the 
U(l)/ charges. However, the quantum numbers in transitions involving the first 
generation are dependent on which of the 3 choices is made. These charges 
cannot be related to a different case by such a linear transformation of the 
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U(l) charges. Clearly the 3 choices correspond to the arbitrary labelling of the 
3 proto-generations according to the U(l)/ charges. There are 3 rather than 
6 choices since changing the sign of the U(l)/ charges simply corresponds to 
interchanging the 2nd and 3rd proto-generations. So, although we will only 
consider the choice made above (a dominant II — cr mass matrix element) in 
this paper, our conclusions about modelling the second and third generations 
will be independent of this choice. 

To begin with we shall concentrate on the mass matrix elements responsible 
for the second and third generation masses. We have chosen the elements in 
the up-type mass matrix which will give the top and charm quark masses. In 
the down-type and lepton mass matrices we must choose the diagonal elements, 
if we want to obtain the order of magnitude degeneracies between the bottom 
quark and tau lepton masses, and between the strange quark and muon masses. 
Therefore, ignoring the elements involving the first generation, we have the 
following order of magnitude elements in the mass matrices: 


Mu 

Mu 

Me 



me \ 
mb ) 

mbV23 \ 
mb ) 
7 

mb 


( 22 ) 

(23) 

(24) 


where ~ ms and mr — mb and we have indicated the likely position of 
the element which leads to the mixing V 23 between second and third generation 
quarks. 

Now we are in a position to choose specific Higgs fields which could lead to 
these types of mass matrices. We can do this by noting that there are relations 
between the quantum numbers of different elements. In particular we define: 

b = QbL - QbR - Qws (25) 


which are the total charges carried by the Higgs fields, other than the WS Higgs 
field, in the transition which leads to the element in the down-type mass matrix 
corresponding to the bottom quark mass. The U(I) charges of the WS Higgs 
field are: 

Qws = Qcn-Qt^ = ( 26 ) 

We define similar quantities for other elements in the down-type and lepton 
mass matrices. In the up-type mass matrix we have a different interaction with 
the WS Higgs field, and so we define the total charges carried by the other Higgs 
fields in, for example, the transition responsible for the charm quark mass by: 

c = QcL-QtR + Qws (27) 
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Here we have used the requirement that the charm quark mass should be due 
to the off-diagonal element corresponding to the transition between cl and Ir, 
rather than the diagonal element corresponding to the transition between cl 
and cr. It is easy to verify the following relations between these charge vectors: 

6 -j- c “t" s = 0 

Qb23 s — b = 0 (29) 

where Qb 23 is defined to be the charges corresponding to the element we ex¬ 
pect to have magnitude mf,V 23 in eq. (p^. This means that these vectors form 
triangles in charge space. Thus we predict that (in units of the supposed un¬ 
suppressed top quark mass) the mass of the bottom, charm or strange quark is 
no less than the product of the other two masses. This is in fact true in nature 
and also for the other case, eq. (p^, with mb, rus and mbV 23 . 

We now try to fit the effective SM Yukawa couplings in terms of as few Higgs 
fields as possible. In the case of only 1 Higgs field W the prediction from eq. ( p^ 
for hs, the smallest Yukawa coupling of the three, is: 

hs ~ hbhc (30) 


To see this note that each Yukawa coupling is of the form 

hx ~ (IH)I”^I 

where from eq. ( p^ 

rib + ric + Us = 0 

and so, since hs < he < hb implies |ns| > \nc\ > |nb|: 

l^sl = |nb| -I- |nc| 


(31) 

(32) 

(33) 


Eq. ( plD shows that eqs. ( |33| ) and (30) are equivalent. The VEV {W) (in units 
oi Mr) is the suppression factor generated each time that the Higgs field W has 
to be applied. Now eq. ( p0| ) is not true in nature (see table |^) and so we must 
use at least 2 Higgs fields. 

Whenever we have relations such as eqs. (p^ and (29) between three matrix 
element charge vectors, we imagine that we can re-use chains of Higgs fields 
(see fig. 1 ) able to exchange the quantum numbers necessary to give masses 
to, say, the b and c quarks to also give mass to the s quark. Naturally, there 
will be a piece of chain common, up to charge conjugation, for b and c not 
used by s. We seek to symbolise this chain structure in figure 2: the connected 
single lines represent logarithms of the suppression factors for the masses of the 
6 , c and s quarks. The lengths of the double lines represent the logarithms of 
the bunches of Higgs field suppression factors common for those two quark 
masses represented by the lines forming the double lines in question. A priori we 
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Hs ~ 4.5 



Figure 2: The magnitudes of the logarithms of the Planck scale Yukawa cou¬ 
plings for b (really t), c and s (really fi) are represented by the total lengths of 
the corresponding single lines. The lengths of the double lines labelled by Hi 
represent — In Hi. 


have the possibility that, even for three linearly dependent transition quantum 
number sets (obeying, say, eq. (|^) like b, c, s), there could be some Higgs fields 
which are used in only one of the three relevant diagrams like figure 1. However 
this is not possible for a set of Higgs fields with linearly independent quantum 
numbers. 

We now assume that just 2 Higgs fields, with linearly independent quantum 
numbers, are needed to give the masses of the b, c and s quarks and denote 
them by W and T. Denoting the quantum number vectors of these fields W 
and T by Qw and Qt respectively, we must be able to write 

b = nbQw+PbQT (34) 

and similar expressions for c and s. Then eq. ( ^^ is satisfied as above and 
similarly 

Pb+Pc+Ps=0 (35) 

So, for each Higgs field, we must have a relation like: 

be I = \pb\ + bsl (36) 

or a similar relation with the quark names permuted. It follows that the sup¬ 
pression factors associated with each Higgs field can be factorised in a similar 
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way to eq. and the 3 Yukawa couplings can be split order of magnitudewise 
in the following way: 


hb ~ niHz (37) 

he 0^ HiHa (38) 

hs n2ll3 (39) 

The IliS are products of Higgs suppression factors, i. e. they are of the form: 

n, = (1T)I”^I(T)IP^I (40) 

Here, and henceforth, we express the Higgs held VEVs (apart from {(l)ws)) in 
units oi Mp. 

Similarly eq. ( p^ requires the splittings: 

hb (41) 

hs ~ n4n6 (42) 

hbV23 ~ HsHg (43) 

In order to motivate our hnal model from the numerology, we choose to 
replace the quark Yukawa couplings hb and hg by the lepton Yukawa couplings 
h^ and h^, because the latter ht our model better. It is anyway our prediction 
that hb ~ hr and hg h^, since the quantum numbers satisfy: 

b = T (44) 

s = jl (45) 

From eqs. (|^)-(|^) or their illustration in hgure 2, we see, by insertion of the 
Planck scale Yukawa couplings from table 1, that: 


In Hi ~ 

-(In /lb + \nhc- In /is) 

(46) 


-(In/iT- + In he - In/i^) 

(47) 


-1.8 

(48) 

lnn2 ~ 

i(ln/i5 + In/is - In/ic) 

(49) 


i(ln/iT- + In/i^ - In he) 

(50) 


-2.8 

(51) 

In Ha ~ 

i(ln/is + In/ic - In hb) 

(52) 


i(ln/i^ + In/ic - \nhr) 

(53) 


-4.5 

(54) 
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Figure 3: The magnitudes of the logarithms of the Planck scale Yukawa cou¬ 
plings for b, s (really /r) and bxV 23 are represented by the total lengths of 
the corresponding single lines. The lengths of the double lines labelled by 11^ 
represent — In 11^. 


Similarly, from eqs. (^)-(|^) or figure we get: 


In 114 

~ -2.1 

(55) 

Inlls 

~ -2.5 

(56) 

In He 

~ -5.3 

(57) 

we have 

a very good order of magnitude 

relation: 

Ha- 

iiin2 

(58) 

' Hi and 

112 as single Higgs field VEVs 

(in units 

(W) 

= Hi 

(59) 

(T) 

= n2 

(60) 


oiMp): 


We also observe that hb — hr is, to a good numerical approximation, split in the 
same way in these 2 cases eqs. (^) and (^). So we make the identifications: 


114 — III 

IIs = 112 
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(61) 

(62) 




This gives us the approximate values of: 


{W) ~ ~ 0.14 (63) 

(T) ~ ~ 0.07 (64) 


It is also numerically suggested to express Hq in terms of the 2 Higgs field VEVs 
as: 


He = (T)2 


(65) 


We can now see that it is consistent to split hg 
have: 


hg = {W){T)^ = {{T)){{W){T))-. 


ns 


in 2 different ways, since we 


((W»((T)2) 



( 66 ) 


So we have shown that it may be possible to fit the masses and mixing angle 
of the second and third generations using the two Higgs fields W and T. In 
particular we make the order of magnitude predictions: 


hg ~ 

{W){Tf 

(67) 

Hq — 

{Wf{T) 

(68) 

hr — hb — 

{W){T) 

(69) 

Vtg-V23 ~ 

(t) 3 (r)2 

hb ~ (W) 

(70) 


Now we have determined all the H^ and so we know the values of \ni\ and \pi\ 
defined in eq. (^). However, we don’t know the signs of rii or pi. Essentially, 
if e.g. m < 0 then \ni\ hermitean conjugated Higgs fields are used for 
the corresponding element in the mass matrix rather than the Higgs fields W 
themselves. By definition we can choose the sign of the charges of W and T so 
that: 

b = Qw + Qt (71) 

Then we can use eq. (|2^), with our reasonable assumption that Qw and Qt are 
linearly independent, to see that the only consistent choice is that 

c = —2Qw + Qt (72) 

s = Qw ~ ‘^Qt (73) 


7 Predictions from a generalised model 

Clearly we can now calculate the charges of the Higgs fields W and T in the 
AGUT model. However, we shall first highlight the predictions that only depend 
on the following general features suggested by the AGUT model and the data 
as discussed in section 4: 
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(a) Left-handed up-type and down-type quarks have the same charges. 

(b) The “diagonal” matrix elements in the up-type mass matrix involve the 
opposite charges to those in the corresponding elements of the down-type 
and lepton matrices. The “diagonal” elements are just meant to be di¬ 
agonal after some appropriate permutation of the imagined left and right 
handed proto-quarks and leptons. 


(c) The top quark mass has zero quantum numbers suppressing it and is “off- 
diagonal” . 


Now the linear relations, eqs. (|28| ) and (M), are in fact consequences of these 
assumptions (a), (b) and (c). Using the phenomenological arguments of the 
previous section, we take as our final assumption: 


(d) The quantum numbers of the Higgs fields W and T should be related to 
those for the mass matrix elements by eqs. 0-®. 

With the above 4 assumptions extracted from the AGUT model and phe¬ 
nomenology, we obtain the Yukawa coupling matrices in the form: 


Hu 


Hd 


He 


1 ) 
f W(Tt)2 r^3 \ 
tU2(yt)4 J 

( W{T^f ? \ 


WT ) 


(74) 

(75) 

(76) 


where the order of magnitudes of the elements are obtained by replacing the 
Higgs fields by their expectation values (in units of Mp)- So only the 2 off- 
diagonal elements in the charged lepton matrix are model dependent. Since we 
now have all the matrix elements for the quarks we can check that we obtain 
the required masses and mixing angle, i.e. that no elements dominate the ones 
we expected to be relevant for calculating the eigenvalues. 

The largest eigenvalue in a 2 x 2 matrix is approximately the largest element 
and the other eigenvalue is approximately the determinant divided by the largest 
element. So we can see that the up-type matrix has the form we expected, 
since no element is larger than 1 and (W)^(T) » ((W)(T)^)((1U)(T)). In the 
down-type matrix we see that the element we wanted to be /it,V 23 has the correct 
magnitude (T)^. We also see that the other off-diagonal element in is greatly 
suppressed and does not give any significant contribution to the eigenvalues. 

As already discussed in section 4, the first generation u, d and e masses 
could be dominated by the extra diagonal elements in the 3 generation Yukawa 
coupling matrices Hp, Hp and Hp- 

It should be noted that when we consider a specific model such as the AGUT 
model, problems can arise. Most obviously the off-diagonal elements in Hp could 
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dominate one or both of the required eigenvalues. In this case the model would 
not be suitable. Secondly, it might turn out that the charges of the Higgs fields 
W and T were not independent and then the elements would be expressible as 
different combinations of these fields. In this case the matrix element would 
correspond to the least suppressed combination and this may be different from 
the elements suggested above. Finally, a similar situation arises when elements 
involving the first generation are considered. Then it may be necessary (as it is 
in fact in the AGUT model) to introduce more Higgs fields, since all the elements 
may not be expressible in terms of W and T alone. If all the Higgs fields are 
not linearly independent then we again have the possibility that there are other 
less suppressed combinations. These last two situations do not necessarily spoil 
the model, but a careful check is required. 


8 Specific choice of Higgs fields within the AGUT 
model 

We are now ready to calculate the charges of the Higgs fields W and T in the 
AGUT model. This is quite simple using eqs. (0) and d?^). We obtain: 

1,- ^ f 1 1 4\ 

Qw = ^{h-<^ = (^0,--,-,--j 

Qr = b-Q^ = (^o,-i,i,-0 

Obviously these charges are linearly independent and so the second potential 
problem, discussed at the end of the previous section, does not arise in this 
model. We can now calculate all the elements involving only the second and 
third generations. The relevant elements of the three Yukawa matrices are given 
by: 


(77) 

(78) 


Hu 

f {W^)^T \ 

V 1 ) 

(79) 

Hd ~ 

( iu(rt)2 J.3 \ 

V W2(Tt)4 WT ) 

(80) 

He ^ 

V (IU1')2t4 WT ) 

(81) 


Now it can clearly be seen that the off-diagonal elements in the charged lepton 
matrix are sufficiently suppressed, so that they do not make significant con¬ 
tributions to the eigenvalues. Thus, as required, the first potential problem 
discussed at the end of the previous section does not arise in this model. So we 
have shown that this model could give realistic fermion masses for the second 


20 




and third generations. We must now complete the model by determining the 
matrix elements involving the first generation. 

We notice that the charges of W and T do not cover the 2 dimensional space 
of charges ^ and yf, since only even yj charges can be constructed with integer 
numbers of these Higgs fields. Therefore, since both W and T have ^ = 0, we 
will need at least 2 more Higgs fields to fully cover the 3 dimensional charge 
space required to break G down to the SMG. We will now choose 2 more Higgs 
fields which, together with W and T, will fully cover this space. 

As a simple proposal, we may introduce a third Higgs field ^ provided with 
charges corresponding to the quantum number differences between the left- 
handed quarks in the first and second generations. Then the Hu and Ho matrix 
elements in the row corresponding to the left-handed first generation quark be¬ 
come the same as those in the second generation row augmented by an extra 
factor The Vus mixing matrix element is then expected to be dominated by 
the Hjj and Ho elements in the first generation row, corresponding to those in 
the second row which dominate the c and s masses. In fact Vus is readily seen to 
be given, in our type of model, by the ratio of these first to second row matrix 
elements and should thus be equal order of magnitudewise to the ^ Higgs field 
suppression factor: 

Vus - iO - 0.22 (82) 

Similarly we expect: 

Vub ^ {OV23 ^ VusVcb (83) 

The quantum numbers for such a ^ field must of course be: 

Qf=Qd^-Qs^ = (^-0,0,0^ - (o,i,o,o) = (^>-^, 0 , 0 ^ (84) 

We must now choose one more Higgs field to fully span the 3 dimensional 
space of charges. Otherwise the first generation would remain massless. In order 
to be consistent with the well-known phenomenological relation: 



which motivated the Fritzsch [p4]| and many subsequent quark mass ansatze, we 
clearly want the first generation “diagonal” matrix elements to be suppressed, by 
a factor (^)^ relative to the Yukawa matrix element dominating the s quark mass. 
With our proposal that the first generation row has a factor ^ more than the 
second generation row, we must thus arrange for the transition matrix element 
from sl to do to have the same factor W(T'^)^^ as that from do to sr. This can 
be achieved—numerically—by proposing a Higgs field S, having a suppression 
factor that is actually equal to unity and that can compensate for the difference 
in quantum numbers between these two hoped for equally suppressed matrix 
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elements. This will lead to 2 different but comparable mechanisms for the down 
quark mass. So we take: 

(5) = 1 (86) 


and the charges of S are given by: 

Qs — \Qsl ~ ~ [Qdh ~ Osr] 


1 


0-g>0,0 ) - ( --,0,0,0 


1 


T,0,0,0 

6 




We note that the SM weak hypercharge vanishes 



(87) 


y = yi + y2 + y3 = 0 


( 88 ) 


for the Higgs fields W, T, ^ and S. This guarantees that the SMG is recovered 
as the diagonal subgroup of the SMGi groups. 

We can now calculate the suppression of all elements in the Yukawa matrices. 
However, we must first note that, since we have used 4 Higgs fields, we cannot 
uniquely resolve the charge differences between left-handed and right-handed 
fermions. There will be some combination of the 4 Higgs field charges which will 
result in vanishing charge differences. We must find the smallest combination 
of the 4 Higgs fields which results in a vanishing set of charges Q = 0. To do 
this we note that all fermion U(l)/ charge differences are quantised as integers. 
However, the 3 Higgs fields W, T and ^ can only give integer U(l)/ charge 
differences which are even. Therefore we must have at least two S Higgs fields 
involved in the combination. Then we can Hnd the simplest combination of the 
other 3 Higgs Helds which, together with the two S fields, give net vanishing 
charge differences (5 = 0. This combination is: 

2Qs - 2 Q 5 - 9Qt + 5Qw = 0 (89) 

Since this involves such large powers of T, there is usually no ambiguity in 
selecting the combination of Higgs fields which suppresses the transition the 
least. So finally we can write out the full Yukawa matrices for each type of 
fermion: 

/ S'1'WtT2(^t)2 (wi')2rc \ 

Hu -- S'1'Wl'T2(^t)3 (^t) 2 y (90) 

V S't(^t)3 I J 

( S'W(Tt)2^2 p^(7^t)2^ rj.3^ \ 

Hd ~ SWiT'^fi W(rt)2 J.3 (9^) 

\ S'W2(Tt)4^ lT2(Tt)4 WT j 

( S'W(rt)2^2 p^(7’t)2(^t)3 (S't)2^r4^t \ 

He ~ SW{T^fi^ W{T^f (92) 

\ S'3W(rt)5^3 WT ) 
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Note that one of the elements in He has been changed after the introduction 
of the Higgs fields ^ and S', but this element is still so suppressed that it has 
practically no relevance for any of the charged lepton masses. 


9 Results 


Now we are able to choose specific values for the 3 VEVs (VH), {T) and (^) and 
calculate the resulting masses and mixing angles. The overall mass scale for 
the ht is set by eqs. (D and (i- In order to hnd the best possible fit we must 
use some function which measures how good a fit is. Since we are expecting 
an order of magnitude fit, this function should depend only on the ratios of the 
fitted masses to the experimentally determined masses. The obvious choice for 
such a function is: 


X 


= E 


In 


/ m 


\ ™exp / _ 


(93) 


where m are the fitted masses and mixing angles and m-exp are the corresponding 
experimental values. The Yukawa matrices are calculated at the fundamental 
scale which we take to be the Planck scale. We use the first order renormalisation 
group equations (RGBs) for the SM to calculate the matrices at lower scales. 
Running masses are calculated in terms of the Yukawa couplings at 1 GeV. 
The only exception is the top quark, where the experimentally measured mass 
is the pole mass and this is what we quote. We present here the result of an 
updated fit, using the following values of the SM gauge coupling constants 
at the electroweak scale and their values extrapolated to the Planck scale via 
the RGBs: 


t/(I): 

gi{Mz) = 0.462 

gi{Mpianck) = 0.614 

(94) 

SU{2) : 

g^iMz) = 0.651 

g 2 {Mpianck) = 0.504 

(95) 

S't/(3) : 

gsiMz) = 1.22 

gi{Mpianck) = 0.491 

(96) 


We cannot simply use the 3 matrices given by eqs. (^)-(^) to calculate the 
masses and mixing angles, since only the order of magnitude of the elements is 
defined. This could result in accidental cancellations if we calculated the eigen¬ 
values and eigenvectors using these values. Therefore we calculate statistically, 
by giving each element a random complex phase and then finding the masses 
and mixing angles. We repeat this several times and calculate the geometrical 
mean for each mass and mixing angle. In fact we also vary the magnitude of 
each element randomly, by multiplying by a factor chosen to be the exponen¬ 
tial of a number picked from a Gaussian distribution with mean value 0 and 
standard deviation 1. 

We then vary the 3 free parameters to find the best fit given by the 
function. We get the lowest value of for the VBVs: 

{W) = 0.179 (97) 
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Table 4: Best fit to conventional experimental data. All masses are running 
masses at 1 GeV except the top quark mass which is the pole mass. 



Fitted 

Experimental 


3.6 MeV 

4 MeV 

rrid 

7.0 MeV 

9 MeV 

me 

0.87 MeV 

0.5 MeV 

me 

1.02 GeV 

1.4 GeV 

ms 

400 MeV 

200 MeV 


88 MeV 

105 MeV 

Mt 

192 GeV 

180 GeV 

mb 

8.3 GeV 

6.3 GeV 

TTlr 

1.27 GeV 

1.78 GeV 

Vus 

0.18 

0.22 

Veb 

0.018 

0.041 

Vub 

0.0039 

0.0035 


(T) = 0.071 (98) 

iO = 0.099 (99) 

The fitted value of (^) is approximately a factor of two smaller than the estimate 
given in eq. (^^. This is mainly because there are contributions to Vus of the 
same order of magnitude from both Hu and Ho ■ The result of the fit is shown 
in table The experimental values are those given in table |^. This fit has a 
value of: 

= 1.87 (100) 

This is equivalent to htting 9 degrees of freedom (9 masses + 3 mixing angles - 

3 Higgs VEVs) to within a factor of exp(-\/1.87/9) ~ 1.58 of the experimental 
value. This is better than would have been expected from an order of magnitude 
fit and should be compared with = 3.7 for the fit with only 7 degrees of 
freedom in . 

We can also fit to different experimental values of the 3 light quark masses 
by using recent results from lattice QCD |^. Light quark masses extracted from 
lattice QCD seem to be consistently lower than the conventional phenomeno¬ 
logical values (H given in table ||. We take the following light quark masses as 
typical lattice values, extrapolated to 1 GeV using the RGBs: 

mu 1.3 MeV (101) 

ma ~ 4.2 MeV (102) 

ms ~ 85 MeV (103) 
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Table 5: Best fit using alternative light quark masses extracted from lattice 
QCD. All masses are running masses at 1 GeV except the top quark mass 
which is the pole mass. 



Fitted 

Experimental 


1.9 MeV 

1.3 MeV 

rud 

3.7 MeV 

4.2 MeV 

me 

0.45 MeV 

0.5 MeV 

me 

0.53 GeV 

1.4 GeV 

ms 

327 MeV 

85 MeV 

m^ 

75 MeV 

105 MeV 

Mt 

192 GeV 

180 GeV 

mb 

6.4 GeV 

6.3 GeV 

m-r 

0.98 GeV 

1.78 GeV 

Vus 

0.15 

0.22 

Veb 

0.033 

0.041 

Vub 

0.0054 

0.0035 


We can now vary the Higgs VEVs to give the best ht to this alternative data. 
The best ht is shown in table ||. The values of the Higgs VEVs are: 


{W) 

(T) 

(0 

and this ht has a larger value of: 


= 0.123 

(104) 

= 0.079 

(105) 

= 0.077 

(106) 

= 3.81 

(107) 


Comparing this ht to the one above, using conventional light quark masses, 
we can see that we have an improvement in the hrst generation masses, since 
the up and down quark masses were lowered towards the electron mass. The ht 
to Vch has also improved. However, the strange quark mass, which we always 
predicted too large, is even worse because it’s experimental mass has also been 
lowered. It may seem that now our assumption, that hg ~ at the fundamental 
scale, is not correct. Of course we must remember that everything should be 
taken order of magnitudewise and we can always ignore one borderline case such 
as this. However, it would be interesting to see if we could have produced a 
model without the order of magnitude degeneracy between the strange quark 
and and muon masses. 
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10 Suggestion for a better model of the alterna¬ 
tive masses 


At first it would appear to be impossible to produce a model without the order 
of magnitude degeneracy hg — h^, while retaining the other order of magnitude 
degeneracies between the leptons and down-type quarks. Indeed, we would cer¬ 
tainly not want to spoil the natural prediction of hr-, which is well-known 
to be quite accurate. The relation between the down and electron masses may 
not appear to be so accurate. However, when we consider the lower value of the 
down quark mass extracted from lattice QCD, along with the hindsight realisa¬ 
tion that we actually predict hd ~ 2he due to the two competing combinations 
of elements giving the lowest eigenvalue, we actually have quite an accurate 
prediction which we would not want to spoil. As we noted in section ^ if one 
eigenvalue involves an off-diagonal term, then so must at least one other. Thus, 
it would appear that we cannot spoil the unwanted relation, hg ~ without 
spoiling another relation. 

However, we can in fact do precisely this. Examining the Yukawa matrices 
given by eqs. and (p^), we see that there is one off-diagonal element which 
has the same order of magnitude in the down-type and charged lepton matrices. 
It turns out that this is not by chance and is actually a consequence of the 
element being in the same position as the unsuppressed element in the up-type 
matrix, which leads to the top quark mass. To see this in general, we use the 
notation of section || for the net charge differences supplied by the Higgs fields 
suppressing each element in the Yukawa matrices. 

If the top mass is to be unsuppressed then by definition we have, for some i 
and j : 

^Quij = QuiL - QujR = QoiL - QujR = 0 (108) 

We now wish to show that, for the same fixed i and j, this implies: 

AQdzj = (109) 


where we know the sign is minus (either sign would lead to the order of magni¬ 
tude equality) because of th e specific example of eqs. (0) and (^ . So now we 
can see that, using eq. (lOS): 


^Qoij = “^QDiL - QojR - QujR ( 110 ) 

-AQsij = -QoiL - QsiL + QujR + Qejr ( 111 ) 

We can see from table ^ that all left-handed fields carry zero U(l)y charge 
and so: 

2(5Dii = -QoiL - QEiL (112) 

since the only non-zero charges are the U(l)i charges, which are equal to the 
weak hypercharges in the SM. An alternative, perhaps more fundamental, reason 


26 



for this equality is that it corresponds to the absence of anomalies. In particular 
the fields DiL (which is the same as UiL) and EiL are the only fields coupling 
to SU(2)i. So eq. ( |112| ) is simply the condition for cancellation of anomalies 
associated with the triangle Feynman diagrams with two external SU(2)i gauge 
bosons and one U(l) gauge boson (with 4 independent choices of the U(l) gauge 
group). 

Similarly the relation: 


- Qdjr - QujR = QujR + Qb.b (113) 


can be seen to be true. Again this is due to anomaly cancellation, though not 
in such a simple manner as the cancellation for the left-handed fields. So we 
now see that eq. (1^) is true. Therefore it may be possible to produce a model, 
where the bottom quark and tau lepton masses come from the element in the 
same off-diagonal position as the unsuppressed element in the up-type matrix. 
This would retain the good relation, /if, ~ /it, as well as the order of magnitude 
degeneracy of the 1st generation masses, but would not enforce the less desirable 
prediction of hs ~ h^. This interesting scenario is currently being investigated 
and will not be further commented on in this paper. 


11 CP violation 

Another prediction, which can be made from a model of the mass matrices, 
is the amount of CP violation due to the CKM matrix. This depends on the 
complex phases in the matrix usually parameterised by a phase J Q. However, 
there are different ways of parameterising the unitary CKM matrix, and so it 
is better to define a parameterisation-independent quantity which is a measure 
of the amount of CP violation. Such a definition is possible and corresponds 
to the areas of the “unitarity triangles” [^. The 3 sides of the triangles are 
defined in the complex plane as Si = VijV*f., where j ^ k are fixed and i labels 
the 3 sides. The condition that the CKM matrix V is unitary determines that: 

Si + S2 + S3 = 0 ( 114 ) 

and so these 3 lines in the complex plane form a triangle. Also, the areas of all 
the different triangles are the same. Therefore we can define the amount of CP 
violation, J, in terms of the area, A, of any of these triangles, e.g.: 

J = 2A = \x 1 y 2 - X2yi\ (115) 

where Xi (yi) are the real (imaginary) components of s^. In our model with the 
conventional values of the light quark masses, we find that the model predicts: 

5.8x10"® (116) 
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In the case where we use the lower lattice values for the light quark masses, we 
have the prediction: 


1.2 X 10"® 


(117) 


These predicted values are both lower than the experimental determination 


m 


J~2.0x 10"®-3.5 X 10 


1-5 


(118) 


The model predictions agree in order of magnitude though, for the case where 
we take the conventional values of the light quark masses, the prediction is 
clearly worse than any of the mass and mixing angle predictions. However, we 
should really even consider this worst case as sufficiently good agreement with 
experiment, since our prediction of the CP violation parameter J is expected to 
be less accurate than that for a typical mass or mixing angle. This is because 
our prediction for J ~ involves 8 Higgs field VEVs, whereas a typical mass 
prediction involves rather of the order of 3 VEVs such as rus ~ WT^. As noted 
after eq. (IOC), with = 1.87, we expect to fit a typical mass within a (one 
“standard deviation”) factor of exp(-\/l.87/9) ~ 1.58 = exp(±0.46). If now the 
logarithm of the uncertainty factor for J is taken to be | times as big, then J is 
expected to b e uncertain by a factor of exp(±1.2) = 3.4^^. S o ev en the worst 
case, eq. (|ll6| ), only deviates from the experimental value, eq. (11^), by a factor 


of log 


2.7xl0~ 


r/1.2 = 1.3 “standard deviations”. Thus we conclude that our CP 


5.8x10“ 

violation predictions agree with experiment, within the accuracy we can expect. 

We will now show that the amount of CP violation in our model can be 
well estimated in terms of the fitted mixing angles. We will first consider the 
simplified case, where the up-type Yukawa matrix is approximately diagonal, 
and so all the quark mixing is due to off-diagonal elements in the down-type 
Yukawa matrix. Then the down-type matrix should be of the form: 



/ hd 

hsVi2 

z 


X 

hs 

hbV23 


V X 

X 

hb 


(119) 


where x denotes element which are considered to be “sufficiently suppressed”. 
This means that the lower off-diagonal components should not be so large as 
to be relevant when calculating masses and mixing angles. Since the matrix 
diagonalised is HdH\, it can easily be seen that these elements will not be 
relevant unless they are considerably larger than the corresponding upper off- 
diagonal elements. This is because they get multiplied by smaller diagonal 
elements, since hd ^ hb. As we will show: 


Vus ^ Wd Vi2 (120) 

Vcb ^ Pt. P 23 (121) 

However, the element 2 ; determines Vub but only contributes to Vtd- If ^ is 
bigger than or of the same order of magnitude as the product htVusVcb, we 
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would expect that we should have: 

Z — hifVub — hfjVtd 


( 122 ) 


but if z were smaller than this value, we would still predict the correct order 
of magnitude of Vtd through the ‘indirect’ mixing between all 3 generations as 
opposed to the ‘direct’ mixing due to z, as we shall now explain. 

Consider how to approximately diagonalise the matrix HdH^. If z is so small 
that there is essentially no ‘direct’ contribution to Vtd, then we can use a 2 stage 
process using unitary matrices of the form: 

U = ( 0 ^ and C = f 0 

V 0 0 1 / Vo 



(123) 


where, in the approximation of small mixing, we have the order of magnitude 
unitary matrices (suppressing all phases for convenience): 




The order these matrices are applied is important. It can be seen that when z 
is small: 

U^V^YdYjvU ~ diagiylvlyl). (125) 

If we had applied the matrices in the opposite order we would not have approx¬ 
imately diagonalised YdY^. Therefore the CKM matrix is given by the product 
VU. With the approximations we have made above, this is: 


/I U2 0 \ 

VcKM — VU ~ I Vi 2 1 V 23 I 

V V12V23 C23 1 J 

So we see that we have the symmetrical relations: 

Vcd — Vus — V12 

Ct, ~ Vcb-V23 

but the mixing between 1st and 3rd generations is given by: 

Vub ^ 0 

Vtd ^ ^12^23 = Vs" 


(126) 


(127) 

(128) 


(129) 

(130) 


This is what we refer to as ‘indirect’ mixing between the 1st and 3rd generations. 

Now we shall consider the case when there is also ‘direct’ mixing between 
the 1st and 3rd generations. This is the case where: 

z ~ hfcV" (131) 
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In this case we can approximately diagonalise HdH\ by using the unitary matrix: 



f 1 

0 

T/direct 

^13 

w= j 

0 

1 

0 


i T/direct 

\ 

0 

1 


(132) 


and then applying the matrices V and U as above. In the approximation of small 
mixing angles, this will only affect the mixing between 1st and 3rd generations. 
We will get: 


Vub ^ 

T/direct 

^13 

(133) 

Vtd 

T/direct . x/indirect 
^13 *^13 

(134) 

So we have found the relation: 



Vtd 

— ^uh ^us^ch 

(135) 


In fact these 3 terms correspond to the 3 sides of one of the “unitarity triangles”. 
The ‘direct’ and ‘indirect’ terms arise from different matrix elements of Hd, 
which have been assigned independent random phases that are averaged over. 
This means that, for given values for the magnitudes of Vub and VusVcb, the 
angle 9 between these 2 sides should be random. We can then estimate the 
amount of CP violation (which is two times the area of this triangle) to be: 

J ~ ^Vu,VcbVub (136) 

where the ^ is the geometric average of sin 0 . 

Actually there is a contribution to the CKM matrix from the up-type Yukawa 
matrix in our model and thus the above considerations are not quite true, in 
the sense that we do not really have then that the phases of two of the sides in 
the unitarity triangle are independent. In fact the combination of the up-type 
12-transition matrix element and the 23-element in the down-type matrix leads 
to a contribution to the CKM matrix, which in turn gives rise to contributions 
to the Vub and the V 12 V 23 sides of the unitarity triangle. By a similar calculation 
as above, and with obvious notation, we have: 


Vtd ^ 

TT-direct . x/dowriTT' 

^13 + ^12 ^23 

(137) 

Vub ^ 

^direct ^ V"iVz3 

(138) 

Vcb ^ 

V 23 

(139) 

Vus ^ 

^down ^up 

(140) 


The phases are such that Vtd, Vub and VusVcb are the three sides of a unitarity 
triangle. It is now clear that there are three independent quantities relevant 
to the triangle: and V^^Vzs- Furthermore, in our model 
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these three quantities are of the same order of magnitude. Therefore there is 
a permutation symmetry among the three sides of the unitarity triangle. With 
such a symmetry it is impossible for the angles to be flatly distributed, since 
each angle must have an expectation value of tt/S in order that their sum be 
TT. However, we expect the distributions to be given by rather smooth functions 
which can naturally be expanded on cos 9, cos 26 etc. (where the angle is called 
9). A contribution of the form of cos(“odd”0) will not influence the average 
of a function of sin0, such as e.g. log sin 0. This is because sine and, thus, 
functions of sine take the same value for an angle and its supplementary angle, 
while cos(“odd”0) changes sign between angle and supplementary angle. So 
only from even cosines such as cos 29, which presumably already comes with a 
rather small coefficient, will we expect any modification of the average of the 
sine of the angle. To the accuracy to which this cos 29 and higher even cosines 
can be ignored we could, for the purpose of evaluating the geometric average for 
sin0, equally well assume the totally flat distribution for 9 used above. Thus 
we do not expect much deviation from the above estimate, eq. (136), of the CP- 
violation strength by the inclusion of the up-type mass matrix contribution. So 
we can still expect a rather good agreement of this estimate with the computer 
calculation. 

We can now calculate our theoretical prediction of CP violation, given the 
fitted values of the mixing angles for the 2 different fits. For the fit to the 
conventional light quark masses we expect: 


J ~ i X 0.18 X 0.018 X 0.0039 ~ 6.3 x 10"® 


(141) 


and for the fit to the alternative light quark masses: 


J ~ i X 0.15 X 0.033 X 0.0054 ~ 1.3 x 10 
2 


-5 


(142) 


We can see that the above predi ctions only deviate by about ten percent 
from the calculated values eqs. (|116| ) and (117). This deviation is due to the 
fact that the up-type Yukawa matrix contributes to the mixing between 1st and 
2nd generations. However, the theoretical prediction still agrees very well and 
is a useful estimate in models where most of the mixing is due to the down-type 
Yukawa matrix and the matrix elements have uncorrelated phases. 


12 Neutrino masses and mixing angles 

We have constructed a successful model of the observed fermion masses and 
mixing angles. It is interesting to investigate what predictions this model would 
give for neutrino masses and mixing angles. We expect the neutrinos to get a 
mass in this model by the usual see-saw mechanism |^, This occurs when 


31 




we treat the SM as a low energy effective theory. Then we can include the 
non-renormalisable interaction: 

Ly = -j^LL^wsHyCiLL^-wsY + h-c. (143) 

IVI p 

with notation as in eq. ( 0 ), where C is the antisymmetric charge conjugation 
matrix and Hy is the effective light neutrino Majorana coupling matrix. After 
electroweak symmetry breaking we obtain mass terms for the neutrinos: 

Cy mass = VZMyCW^ + h.c. (144) 


where the Majorana mass matrix: 


My = Hy 


{^WS? 

2Mf 


(145) 


is necessarily symmetric. Such terms arise in our model in the same way as the 
Yukawa terms for the other fermions. The only difference is that the Feynman 
diagrams corresponding to fig. 0 involve 2 WS Higgs fields and the transitions 
are between and 

The Majorana mass matrix My is, of course, a symmetric matrix but other¬ 
wise, in models of our type with approximately conserved chiral C/(l) charges, 
the matrix elements are generally of different orders of magnitude. As pointed 
out in 1^ , there are two qualitatively different types of eigenstate that can arise 
when diagonalising a symmetric mass matrix with such a hierarchical structure. 
In the first case, a neutrino can dominantly combine with its own antineutrino 
to form a Majorana particle with small mixing angles. The second case occurs 
when a neutrino combines dominantly with an antineutrino, which is not the 
CP conjugate state, to form a 2-component massive neutrino. For example the 
bare r-neutrino might combine with the bare /i-antineutrino. Such states au¬ 
tomatically occur in pairs, with order of magnitudewise degenerate masses and 
maximal mixing. In the example given, the other member of the pair would be 
formed by combining the bare /i-neutrino with the bare r-antineutrino state. 
In particular, this happens when the matrix element with the largest order of 
magnitude is off-diagonal. 

We define the lepton mixing matrix analogously to the CKM matrix. So we 
find unitary matrices Ue and Uy such that: 

uI^MeMIUe = diag{mg,TO^,m^} (146) 

UlMyMlUy = (147) 

Then the lepton mixing matrix is defined by: 

u = uIUe (148) 
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Table 6: Constraints on neutrino mass squared differences Am^ and mixing 
angles sin^ 20 from fits |l^ to the solar and atmospheric neutrino data. 


Experiment 

Mixing 

Am^(eV^) 

sin^ 29 

Solar (Vacuum) 

^/i,T 


> 0.7 

Solar (MSW) 


4 X 10-6 - 10-6 

3 X 10-6 - 10-2 

Atmospheric 

- k'r 

5x 10-6-3X 10-2 

> 0.65 


By this mechanism we expect small neutrino masses, which may just be 
large enough to be observable (if the electron neutrino is one member of an 
almost degenerate mass neutrino pair and hence has a large mixing angle) as 
solar neutrino vacuum oscillations. Our choice of the Planck scale as the fun¬ 
damental scale, Mp = Mpianck, would require the element responsible for the 
electron neutrino mass in the Major ana coupling matrix H^, to be essentially 
unsuppressed (i.e. of order unity). In this case ~ k ~ ^ ^ 10“® eV, 

as appropriate for “just-so” (energy-dependent) vacuum oscillations |]^ There 
is also some experimental evidence for an oscillation between muon and tau 
neutrinos. This would explain the reduction in the number of muon neutri¬ 
nos compared to electron neutrinos, observed on the ground after cosmic ray 
interactions in the atmosphere. 

There are 2 alternative mechanisms for solar neutrino oscillations. There is 
the obvious possibility of vacuum oscillations, i.e. oscillations between the sun 
and the earth |l^ . In order to obtain an oscillation probability which depends 
on the energy of the solar neutrino, as appears necessary when comparing the 
data to the conventional solar model calculations 0,0, the neutrinos reaching 
the earth must have undergone approximately one oscillation. The correspond¬ 
ing “just-so” neutrino mass squared difference is Am? ~ 10“^^ — 10“^° eV^. 
The well-known alternative to the large mixing vacuum oscillation solution is 
the MSW effect |]^. This is essentially an enhancement of small mixing os¬ 
cillations, due to electron neutrinos interacting with electrons within the sun. 
Current solar and atmospheric neutrino experiments [0 constrain the difference 
of neutrino masses squared and mixing angles as shown in table 

However, we note that results from recent solar model calculations |0] give 
a predicted neutrino flux which varies by more than a factor of two. If 
the flux is taken to be an adjustable parameter within this range, it seems 
possible to get an acceptable energy independent fit to all the solar neutrino 
data 10] . By an energy independent solution we mean that is sufficiently 
large that the Ve oscillation occurs many times between the sun and the 

earth, and what is observed is the average probability of oscillation, which does 
not depend on the energy of the neutrino. Such an energy independent large 
mixing vacuum oscillation solution to the solar neutrino problem corresponds 
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to a mass squared difference 10“^ eV^ > Am^ > 10“^° eV^. The upper limit is 
provided by reactor oscillation experiments ||^ . 

We can calculate the neutrino mass matrix in our model, if we assume that 
there are no new Higgs fields which were not involved in the other fermion mass 
matrices. In terms of the Higgs fields already introduced we find: 

/ (S't)2(Wl')2T4(^t)4 (5't)2(p^At)2y4^t \ 

Hu- (5't)2(w1')2r4^t W{T^f (ITl')2r 

\ (Wl')2T(^t)3 {W'^fT 5'2(TTt)2(rt)2(^t)2 j 

(149) 

where, as usual, we assume that all fundamental Yukawa couplings are of order 
1. Clearly all the elements of H^, are suppressed. The largest element is off- 
diagonal and of order (IT)^(T). As emphasised above, the matrix is symmetric 
and so we obtain a pair of almost degenerate eigenvalues. We find, using the 
fitted values eqs. ( p7[]^ , the following eigenvalues for Hu'. 

hu^ ^ hu^ - {W)^{T) - 2.3 X 10-3 (150) 

hu, - {W)^{T)\Cl^ - 7.8 X 10-“ (151) 


Since a pair of off-diagonal elements dominate the matrix, there is almost 
maximal mixing (sin^ 29 = 1) between and with very small mixing 
(sin^ 29 — ~ 10-®) of Vg. This is not suitable for observable solar neutrino 

oscillations. However, we do have the correct mixing structure for atmospheric 
neutrino oscillations. The problem here is that this would require the difference 
in masses squared of and Vr to be: 

Aml^ = \ml^ -mlj = Aml^p - IQ-^ eV^ (152) 


In our model the main contribution to the mass difference, between the quasi¬ 
degenerate mass eigenstates and comes from the largest diagonal element 
in the 2 x 2 s ubmatrix of H^, containing the dominant off-diagonal elements. 
From eq. (1^), we see that this element is and hence: 


~ (r)(^)2 ~ 7 X 10-4 

hu.. 


(153) 


Unfortunately our scale is totally wrong, since with the natural choice of the 
Planck scale as the fundamental scale we predict: 

TO. ~ - 7 X 10-® eV, ruu, - 2 x IQ-^® eV (154) 

2 Ad F 

and 

Ato^ = 2^^to 2 ~ 7 X 10-2° eV2 (155) 

mu, " 


34 







Hence if this model is to generate observable mixing it will be necessary to 
introduce a new mass scale {Mp ^ 10^^ GeV), which although both arbitrary 
and aesthetically unappealing is usually necessary in models of neutrino mass. 
Altering our interpretation of as the Planck scale would spoil the AGUT 
predictions of the gauge coupling constants in j^, and we are reluctant to do 
so. However, it should be noted that an SMG 0 U{1)^ model with the same 
U(l) charges as our AGUT model would be anomaly free. Such a model would 
have the same fermion mass matrix structure but would not suffer from the 
same objection to a lower Mp, as it would have no prediction for the values 
of the gauge coupling constants. This SMG® C/(l)^ model does not appeal so 
strongly to us, precisely because it does not predict values for the gauge coupling 
constants and the choice of the U(l) charges seems rather arbitrary [Q. 

We may introduce an ad hoc new mass scale into the AGUT model and keep 
Mp = Mpianck, by adding another Higgs particle (A say) with an arbitrarily 
chosen VEV, which is a triplet under SU(2) in the SM. Then including the 
interaction: 

= Llt.AH^ClI + h.c. (156) 

where r are the Pauli spin matrices, leads to the neutrino mass matrix: 

M, = H,{A°) (157) 

Here A° is the neutral component of the Higgs triplet. The model can now 
give suitable masses and mixing for the atmospheric neutrino problem, by the 
assignment of the following quantum numbers, corresponding to the exchange 
of 2 WS Higgs helds, to A: 

Qa = -2Qws = ^ 0 , 3 ’ “ 2 ^ ( 158 ) 


It is taken to be a triplet under SU(2) in the SMG 3 group, but otherwise to 
belong to singlet or fundamental non-Abelian representations of the three SMGi 
groups. With this A we have the same prediction, eq. (1^), for 77,^, but with 
the choice of a new mass scale (A°) ~ 1300 eV we now have: 


Am^ ~ 10-2 eV^ 
ruv ~ ~ 3 eV, ~ 10“^ eV 


(159) 

(160) 


So, as well as explaining the atmospheric neutrino problem, the neutrinos 
and Vt would also be hot dark matter candidates. There is also some evidence 
for X/g neutrino oscillations from the LSND collaboration |^; however, 

as we shall see, we would not predict any observable mixing there. The lepton 
mixing matrix U for our model is given by: 



/ 

1 

(0^ 

{T?{0 \ 

/ 1 

10-3 

4 X 10-3 

u ~ 


(«>" 

1 

72 


- -7 X10-; 

0.71 

-0.71 


, 

(«>" 

1 

J_ i 

V 7 X 10-4 

0.71 

0.71 


\ 

V2 

V2 

V2 / 

\ 




( 161 ) 
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where we have used the fitted VEVs of eqs. (^-99) and ignored CP violating 


phases. Here Uia denotes the mixing between flavour eigenstate I and mass 
eigenstate a. This leads to the bound: 


Pl7 


< 




— / ^ \Uu^aUv^l3Uu^aUv^j3\ 

a,/3 


~ 2 X 10 


-6 


(162) 


on the probability for ^ Fg oscillations. This probability is too small to be 
observed by current experiments and, in particular, is incompatible with the 
LSND observation for which: 


Pl7 


~ 10 


-3 


- 1 eV 


(163) 


The above choice of quantum numbers for A is not the only possibility 
and, by choosing appropriate charges, it is in fact possible to explain the solar 
neutrino problem. For example the choice: 


Qa ( 2’ 2’®’3 


(164) 


for the Abelian charges, and singlet or fundamental non-Abelian representations 
under the SMGi groups, leads to the mass matrix: 

/ S'(rt)3lT(^t)4 SiT^fWS} \ 

- (A°) SiT'^fWS,^ (165) 

\ S'r3(lFl')2^t ) 

In this case the largest element is diagonal and corresponds to the mass 

of Vt, while the masses of the lighter and t'g states are quasi-degenerate 
corresponding to the off-diagonal elements {M„)i 2 = For (A°) ~ 3 eV 

we find that 


Am% ~ 7 X 10”“ eV^ 


~ TO. ~ 2 X 10"^ eV, 


~ IQ-i eV 


(166) 
(167) 
is suitable for 


which, since we have almost maximal mixing between Vg and 
the “just-so” vacuum oscillation solution to the solar neutrino problem given in 
table ^ 

As we noted earlier, the allowed range of values for sin^20 and Arn^ from 
table ^ is for energy dependent solutions of the solar neutrino problem. However 
it is also possible to get energy independent solutions with maximal mixing and 
ArUg^ > 10“^° eV^. In part icular this means that we can increase (A°) in the 
above mass matrix, eq. (165) to give: 

(A°) ~ 100 eV (168) 

Ato^ ~ 8 X 10"® eV^ (169) 


~ TO. ~ 6 X 10 eV, 


TO. 


~ 3 eV 


(170) 
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Then we have an energy independent solution to the solar neutrino problem, 
with a sufficiently heavy Vr to provide a candidate for hot dark matter. 

We conclude that it is possible to obtain observable neutrino oscillations 
in our model, but only at the expense of introducing a new mass scale. It is 
quite natural to obtain a pair of neutrinos with quasi-degenerate masses and a 
maximum mixing angle of 0 = 7r/4. Otherwise all mixing angles are expected 
to be small. This is true for a rather general class of models which generate 
the hierarchical structure of the effective low energy Majorana neutrino mass 
matrix by a similar mechanism to ours, using new approximately conserved 
chiral (gauge) charges. In particular we do not expect to obtain signihcant 
mixing between all three neutrinos. 


13 Conclusions 


We have described the development of our proposed system of Higgs fields, 
which generates a realistic quark-lepton mass spectrum and breaks the anti¬ 
grand unihcation (AGUT) gauge group SMG^ down to the SM gauge 

group SMG. The most important feature used in this development is that, in our 
AGUT model, the diagonal elements are suppressed to the same degree in each 
of the three charged fermion mass matrices Mu, Mo and Me- By using this 
property and consideration of possible values for the Abelian gauge quantum 
numbers, we were led to a scheme in which the Planck scale Yukawa couplings 
for the two heaviest generations were fitted in terms of two suppression factors 
W and T. A couple of numerical coincidences supported this picture and we 
ended up with the suggested formulae, eqs. 
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These expressions could arise from a general quantum number system, shar¬ 
ing a few properties with our AGUT model. The most important property is 
the rule that the same quantum number differences suppress the diagonal ele¬ 
ments for a given generation in each of the three mass matrices. Also the top 
quark mass is assumed to be unsuppressed by the quantum numbers. Finally it 
has to be checked that the other matrix elements, assumed to be small in the 
phenomenologically suggested Yukawa coupling matrices in terms of W and T, 
do not come out to be too large. It turned out that our AGUT model could 
indeed realize the general scheme, in terms of two Higgs fields vacuum expecta¬ 
tion values (IT) and (T) expressed in units of the fundamental (Planck) mass 
scale. 

Then, more trivially, the model was extended to include the first generation; 
with all three first generation particles predicted to have order of magnitudewise 
equal masses. This was done by introducing two more Higgs helds, ^ and S, 
the latter of which however gives no suppression, i.e. we took {S) = 1. That 
really means we could have left out the Higgs held S and replaced the AGUT 
gauge group SMG^ ® U{l)f by the subgroup SMG 12 (S) SMG 3 0 C/(l), which 
survives the spontaneous breakdown due to S. We managed to naturally include 
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the well-known order of magnitude relation 14s — followed 

that the mass expected for the d quark should be larger, by a factor of order 2, 
than its first order approximation of being degenerate with the u quark and the 
electron at the Planck scale. 

We were thereby led to the choice of quantum numbers given in eqs. ( [T^ , 
[^ , 1^ , ^7|) for the Higgs fields suppressing the masses, and in eq. (|^) for the 
Weinberg-Salam Higgs field of our model. So finally the charged fermion mass 
matrices are given, order of magnitudewise, by the Yukawa matrices Hu, Hd 
and He oi eqs. 

A fit to the quark-lepton masses and mixing angles was made, using a com¬ 
puter program to insert and average over random complex numbers of order 
unity independently multiplying each of the matrix elements of Hu, He and 
He- The worst deviations from experiment are the strange quark mass, pre¬ 
dicted to be around 400 MeV or twice its conventional value, and a value of the 
mixing matrix element 14h — 0.018 around half its experimental value. Since 
our model only pretends to make order of magnitude predictions, even these 
worst cases should be considered as agreement with experiment. We also made 
a fit using a preliminary lattice QCD determination of the light quark masses. 
Since these lattice values are significantly smaller than the conventional values, 
and in particular we take the s quark mass to be 85 MeV, the second fit is 
markedly worse; it doubled the average square logarithmic deviation between 
our predictions and the data. This result is a rather direct consequence of the 
simple prediction of our model that the order of magnitude of the muon mass 
and the strange quark mass should be the same at the Planck scale. 

We have also estimated the amount of CP violation in our model. Since CP 
violation arises as a product of rather many of our suppression factors, we expect 
our prediction to be more uncertain than that for a typical mass or mixing angle. 
In fact we predict the Jarlskog invariant J measuring CP violation to be a factor 
of 3 or 4 below the experimental value. With our expected uncertainty in the 
exponent, this corresponds to approximately one “standard deviation”. 

Finally we considered the extension of our model to the neutrino sector. The 
natural see-saw mass scale is the Planck mass in our model, for which the only 
possible observable effect would be “just-so” solar neutrino vacuum oscillations 
provided the electron neutrino mixes strongly with or Vt- The symmetri¬ 
cal Majorana neutrino mass matrix, being hierarchical, can naturally give rise 
to pairs of quasi-degenerate neutrino mass eigenstates with maximum mi xing . 
In our model, the effective light neutrino Majorana coupling matrix eq. ( |l49[) 
has this property, but with maximal mixing between and Vt, while Ve be¬ 
comes much lighter than and mixes very weakly with the other neutrinos. So 
we predict no observable neutrino oscillations, unless we modify our model and 
introduce another mass scale into the theory. We have discussed some ways of 
introducing a new mass scale into our model and given examples with observable 
phenomenological implications. These all give neutrino oscillation phenomenol- 
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ogy corresponding to a pair of neutrinos with quasi-degenerate masses and a 
mixing angle 0 = ^. 

Apart from the neutrino sector, the AGUT model gives a good overall order 
of magnitude fit to the fermion masses and mixing angles. 
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